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Abstract. We discuss the geometric phases and flux densities for the metastable states of hydrogen with 
principal quantum number n = 2 being subjected to adiabatically varying external electric and magnetic 
fields. Convenient representations of the flux densities as complex integrals are derived. Both, parity con- 
serving (PC) and parity violating (PV) flux densities and phases are identified. General expressions for the 
fiux densities following from rotational invariance are derived. Specific cases of external fields are discussed. 
In a pure magnetic field the phases are given by the geometry of the path in magnetic field space. But for 
electric fields in presence of a constant magnetic field and for electric plus magnetic fields the geometric 
phases carry information on the atomic parameters, in particular, on the PV atomic interaction. We show 
that for our metastable states also the decay rates can be influenced by the geometric phases and we give a 
concrete example for this effect. Finally we emphasise that the general relations derived here for geometric 
phases and flux densities are also valid for other atomic systems having stable or metastable states, for 
instance, for He with n — 2. Thus, a measurement of geometric phases may give important experimental 
information on the mass matrix and the electric and magnetic dipole matrices for such systems. This could 
be used as a check of corresponding theoretical calculations of wave functions and matrix elements. 
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PACS. 03.65.Vf Phases: geometric; dynamic or topological - ll.30.Er Charge conjugation, parity, time 
reversal, and other discrete symmetries - 31.70.Hq Time-dependent phenomena: excitation and relaxation 
processes, and reaction rates - 32.80.Ys Weak-interaction effects in atoms 



1 Introduction 

In this paper we study properties of geometric phases and 
geometric flux densities for metastable hydrogen atoms in 
external electric and magnetic fields. Geometric phases in 
quantum mechanics were introduced in [1] and have been 
studied extensively since then; see for instance [2-4] and 
references therein. For a discussion of geometric phases 
for systems described by a non-hermitian Hamiltonian see 
[5-10] and references therein. In our group the adiabatic 
theorem and geometric phases for metastable states were 
studied in [11,12]. Both, parity conserving (PC), and par- 
ity violating (PV) geometric phases were identified. One 
aim is to apply the theory developed in this way to the 
measurement of parity violation in light atoms like hydro- 
gen with the longitudinal spin echo technique; see [13-15]. 
But, clearly, a measurement of geometric phases is very 
interesting by itself since these phases represent a deep 
quantum-mechanical phenomenon. For metastable states 
these phases are complex and, therefore, geometry also in- 
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fluences the decay rates of these states, as we shall demon- 
strate explicitly below. 

In the present paper we are primarily interested in the 
structure of PC and PV geometric phases and flux densi- 
ties, that is, what one can say on general grounds about 
their dependences on the external electric and magnetic 
fields. Since a measurement of PV geometric phases is 
one possibility to study atomic parity violation (APV) 
we briefly refer to recent work discussing the present sta- 
tus of this field. Standard reviews of APV can be found 
in [16,17]. A very recent survey of the past, present, and 
prospects of APV is given in [18]. Experimental results for 
the heavy atoms Cs [19-21], Bi [22], Tl [23,24], Pb [25], 
and Yb [26] have been published. See also the review in 
[27]. A large effort is being undertaken to measure APV in 
Ra+ [28,29], and plans for the future FAIR facility at GSI, 
Darmstadt, include a program of APV studies for highly 
charged ions [30-34] . The situation for APV in the lightest 
atoms, H and D, is nicely summarised in [35,36]. In these 
latter papers it is also stressed that from the theory point 
of view H and D are the ideal candidates to study APV. 

Our paper is organised as follows. In Section 2 we in- 
troduce the atomic systems which we want to study. In 
Section 3 we define the geometric phases and flux densi- 
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tics and derive useful representations for them in terms 
of complex integrals. Section 4 is devoted to a study of 
the structure of these phases and flux densities following 
from rotational invariance. In Section 5 we discuss specific 
cases. Section 6 presents our conclusions. In Appendix A 
we explain the notations used throughout our work and 
provide many useful formulae as well as essential numer- 
ical quantities. In addition, we give the non-zero parts of 
the mass matrix for the n ~ 2 states of hydrogen. In the 
Appendices B, C and D we present detailed proofs for 
the relations derived in Sections 3, 4, and 5, respectively. 
In Appendix E (online only) we give explicit formulae for 
various matrices used in our paper. If not stated otherwise 
we use natural units with h = c = 1 . 



we can apply the standard Wigner-Weisskopf method [39, 
40] . The derivation of this method and its limitations are 
discussed in many textbooks and articles, see for instance 
[41-47]. A derivation starting from quantum field theory 
can be found in [48,49]. Let us note that for more com- 
plex situations than discussed in the present paper, for 
instance if radiative transitions are induced between the 
n = 2 states by an external field, we would have to use 
other methods, master equations, the optical Bloch equa- 
tion etc.; see [42]. 

Thus, for the situations we are considering the basic 
theoretical tool is the effective Schrodinger equation de- 
scribing the evolution of the undecayed states with state 
vector |f) at time f, in the Wigner-Weisskopf approxima- 
tion, 



2 Metastable hydrogen states in external 
fields 

We are interested in the states of hydrogen with principal 
quantum number n = 2. Their energy levels in vacuum 
are shown in Figure 1. The lifetimes r of the 2S and 2P 
states of hydrogen in vacuum are ts = Fg^ ~ 0.1216s 
and Tp = Fp^ = 1.596 x 10~^s, respectively; see [37,38]. 
Here Fs^p are the decay rates. We have 16 states with 
n — 2 for which we use a numbering scheme a = 1, . . . , 16 
as explained in Appendix A, Table 2. 

In this paper we shall consider n = 2 hydrogen atoms 
at rest subjected to slowly varying electric and magnetic 
fields. In vacuum the 2S states of hydrogen are metastable 
and decay by two-photon emission to the ground state. 
The 2P states decay to the ground state by one-photon 
emission. Energetically allowed radiative decays from one 
n = 2 state to another one are completely negligible. This 
remains true when we consider the n = 2 states in an 
external, slowly varying, electromagnetic field in the adia- 
batic limit. There, by definition, the variation of the exter- 
nal fields has to be slow enough such that no transitions 
between the n = 2 levels are induced. In this situation 
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Fig. 1. Energy levels of the hydrogen states with principal 
quantum number n = 2 in vacuum. The numerical values of 
the fine structure splitting A, the Lamb shift L and the ground 
state hyperfine splitting energy A are given in Table 1 of Ap- 
pendix A. 



i|lf)=^(5(i),5(t))li) 



(1) 



Here 



^{£{t), Bit)) =Aq-D- e{t) - fi ■ B{t) (2) 



with ^y£_Q the mass matrix for the n = 2 states in vacuum, 
D and /x the electric and magnetic dipole operators, re- 
spectively, in the n = 2 subspace, and £ and B the electric 
and magnetic fields. We are interested in parity conserv- 
ing (PC) and parity violating (PV) geometric phases. The 
mass matrix ^ q will therefore be split into the PC part 
^ Q and the PV part 6 ^py . 



Jf-n + 5^ 



PV 



(3) 



In the standard model of particle physics (SM) 5 ^py is 
determined by Z-boson exchange between the electrons of 
the hull and the quarks in the nucleus. Here, as in [12], we 
split off a (very small) numerical factor 5 from the PV part 
of the mass matrix characterising the intrinsic strength of 
the PV terms. In Appendix A we give explicitly 5. The 
matrices ^ q. D and fj. are discussed further in Appendix 
A, and their explicit forms used for numerical purposes 
are given in Appendix E; see Tables 3, 4, and 5. 

In the presence of electric fields the metastable 2S 
states will get a 2P admixture making them decay faster, 
see Figure 1 of [12]. We are interested in the situation 
where the lifetime of the metastable states is still at least 
a factor of 5 larger than that of the other states. As shown 
in (27) of [12] this limits us to electric fields 



\£\ <250V/cm. 



(4) 



In [11, 12] the adiabatic theorem and geometric phases 
for metastable states were studied and in the present work 
we shall apply and extend the results obtained there. The 
mass matrix ^^ in (1) depends on the slowly varying pa- 
rameters £ and 3. Thus, we have a six-dimensional pa- 
rameter space. Geometric phases are connected with the 
trajectories followed by the field strengths as function of 
time in this space. 
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In the following we shall, for general discussions, de- 
note £ and B collectively as parameters K ^ 



(5) 



Indices i,j £ {1,2,3} will be normal space indices, for 
instance £i, Bi etc., or refer to any three particular com- 
ponents of K . Indices a, b shall refer to the components 
Ka, a € {1, . . . , 6}. The mass matrix (2) shall be consid- 
ered as function of the six parameters K = Ka 





(Ki\ 




K2 


£\ 


Ks 


b)~ 


Ki 




K5 




\kJ 



^{£,lS) = .JaK) 



(6) 



We shall assume that we work in a region of parameter 
space {K space) where ^{K) can be diagonalised. There 
are then 16 linearly independent right and left eigenvectors 



(7) 



\a,K)^E^{K)\a,K) 

(a = l,...,16). 

These eigenvectors satisfy 

{a,K\(3,K) = Sap . 

As normalisation condition we choose 

{a,K\a,K) = 1 

(no summation over a) . 

The complex energies are 



Ea{K) ^ EaR{K) - -r^iK) 



(8) 



(9) 



(10) 



with EaR the real part of the energy and FaiK) = 

(to:{K)) the decay rate, that is, the inverse lifetime of 
the state \a,K). In the following we shall suppose 



\Ea{K)-Ef3{K)\>c>0 



(11) 



for all a j^ /3 where c is a constant. Exceptions where (11) 
is not required to hold will be clearly indicated. 

Below we shall make extensive use of the quasi projec- 
tors defined as 



TPaiK) = \a,K){a,K 
These satisfy 

¥a{K)I>f){K) : 



¥aiK) for a = /3 
for a 7^ /3 



(12) 

(13) 
(14) 



but, in general, the TPa{K) are non-hermitian matrices. 
Furthermore, we shall need the resolvent (C — ^(K)) 
where ^ is arbitrary complex. With the help of the quasi 
projectors we get 

(C - MK)) "" - E (^ - ^" (^)) ~"P" (^) ' 

a 

(n = 0,l,2,...). (15) 

Now we come back to the effective Schrodinger equa- 
tion (1) which reads, replacing {£, B) by K, 



i^^\t) = MK{t))\t) . 

The state vector is expanded as 

16 

|i) = E^"WI"'^W) • 



(16) 



(17) 



We always suppose slow enough variation of the parameter 
vector K{t). As shown in Section 3 of [12] we get then the 
solution of (16) for the metastable states as follows. 
We consider an initial metastable state at time t = 



|t = 0) = EV'a(0)|a,/^(0)) 



(18) 



ae/ 



where a only runs over the index set of the metastable 
states, 

/ = {9, 10,11,12} , (19) 

in our numbering scheme. Then we have, for i > 0, 

\t)^Y.^a{t)\a,K{t)) (20) 

where 



Ipaii) = exp [ - iipa{t) + i7a(0] V'a(O) i 

Mi)= I dt'Ea{K{t')) , 
Jo 

7a(t)= fdt'{a^t')\i-^\a,K{t')) 

<J 



dt' 



(21) 
(22) 

(23) 



The quantities ipa(t) and 7a (i) are the familiar dynamic 
and geometric phases, respectively. For metastable states 
both will in general have real and imaginary parts. 

Here and in the following the labels a = 9, 10, 11, and 
12 correspond to the states \a,K) = \2Si/2,F,F3,£,B). 
These originate from the states 1251/2, F, F3) with (F, F3) = 
(1,1), (1,0), (1,-1), and (0,0), respectively, through the 
mixing with the 2P states according to the PV, the £, and 
B terms in the mass matrix (2). This numbering has to be 
carefully defined, see Appendix A, since we have to follow 
the states in their adiabatic motion along trajectories in 
parameter space. As explained in Appendix A {F, F3) are 
then only labels of the states, no longer the total angu- 
lar momentum quantum numbers. Thus, in order to avoid 
confusion, we shall stick to the labels a for our states in 
the following. 

Below we shall study in detail the geometric phases for 
metastable states for the case that K{t) makes a closed 
loop in parameter space. 
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3 Geometric phases and flux densities 

In this section wc shall discuss general relations and prop- 
erties for geometric phases and the corresponding flux den- 
sities defined below. These relations hold for any system 
with time evolution described by an effective Schrodinger 
equation 



Here we use 

(d(^|) 1/3, K) + (M?|d|/3, K)=0 



(31) 



i§-^\t)^MKit)m 



(24) 



which follows from (8). Note that in (30) a is the index of 
a metastable state, a G /, but in the sum over /? all states 
with f3 y^ a have to be included. 

As a further relation following directly from (29) we 
get the generalised divergence condition 



with NxN matrices ^^^{K), and having metastable states. 
The parameter vector K can have any number of compo- 
nents and the dependence of ^^ on K need not be linear 
as for ^{S, B) = ^{K) in Section 2. 

We consider now the system over a time interval 
< t < T where the parameter vector K{t) runs over a 
closed curve C 



which implies 



idd{a,K\d\a,K)=^0 (32) 



d 
dK'a 



YaMK) 



d 
dKi 



d 



Yo,,UK) + ^r^Y^^ab{K) = 
aKc 



(33) 



C : t-^ K{t), te [0,T], K{T) = A'(0) 



(25) 



The geometric phases (23) acquired by the metastable 
states are then 

7„(C)^7a(r)= / dt'{aJ<rt')\i^\a,Kit')) 



dt' 



= ia,K\id\a,K) , 
Jc 



(26) 



a € /, where / is the index set of the metastable states. 
For our concrete hydrogen case / is given in (19). Here and 
in the following we use the exterior derivative calculus; see 
for instance [50]. Let J-" be a surface with boundary C, 



Let us now derive further representations and properties 
for the geometric flux densities. From (7) and (8) we get 

for /3 7^ 7 

{PJ<\MK)\l,K)^0. (34) 

Taking the exterior derivative in (34) gives 

[EpiK)~E^iK)]{^\d\j,K) 

+ (/Mv|(d^(i^))|7,A-) = 0. (35) 

Since we suppose (11) to hold for a\\ K E J^ wc get for 

/3 7^7 



dJ' = C , 



(27) 



{/3,K\d\j,K) 



{(i.K\{dMK))\l,K) 



and suppose that ^^(K) can be diagonalised for all K E 
T, and (11) holds for the eigenvalues. We get then 



Ep{K)-E,{K) 



(36) 



Inserting this in (30) gives 



7„(C)=i/ {a,K\d\a,K) 
= i / d(^\d\a,K) 

= / Y^^ab{K)dKahdKt 



Yc,,ab{K)^^Y.^E^(^)~Ep{K)] 



P^a 



x(.7^l^#^l,,^)(,7^|Mm,.,^.) 



(28) 



dKa 



dKu 



(a <-> 



Here we define the geometric flux densities Ya^ab{K)^ the 
analogues for the metastable states of the quantities V of 

[1], by 

Ya^abdKaAdKb = \d{aJ<:\d\a,K) , 

Ya,,ab{K)+Y^.ba{K)=0 ■ (29) 

From (29) we get easily 

r«,ab(/i ) dKa A dKb = +i(d(M?|) A (d|a, K)) 

= -i ^ (M?|d|/3, K) A (/Mi'ldla, K) . 



l-Y,[E^{K)~Ep{K)] 



P^a 



X Tr 



[Pa(if) 



dMK)^ ^j^^dMK) 



dKa 



dKb 

(a ^ 



(37) 



/3#« 



(30) 



where we use the quasi projectors (12). 

We shall now derive an integral representation for 
Ya,ab{K). Consider the complex C, plane, see Figure 2, 
where we mark schematically the position of the energy 
eigenvalues Ea{K) and Ej3{K), (3 ^ a. Since we suppose 
(11) to hold we can choose a closed curve Sa which en- 
circles only Ea{K) but where all Ep{K) with (3 ^ a are 
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X 



Fig. 2. The complex (^ plane with the (schematic) location of 
the energy eigenvalues Ec,{K) and Ep{K) with /3 7^ a. The 
curve Sa encircles only Ea{K). 



outside. The geometric flux densities (29), (37) are then 
given as a complex integral 



I Zni 



dCTr 



s„ 



1 



d^iK) 



iC-MK)) dKa 
1 dMK)- 



{C,~MK)f dKb 



(38) 



fi^a 



X Tr 



Pa(/0 



dKadKb 



Tp{K) 



dK, 






OK, '' ' dKadKb 
3 



Y, [Eo^iK] - EpiK)] 



p^a 



xTr[-2P„(X) 






d^{K) 



3K, "f""^'' 3K. 



P,(A-,^#^P„(A-)^#^F„(A-)"=*'*'' 



dKa 



dKb 






dK, 
d^{K) 



dKa "' ' dKb "' ' dK, 



- J2 [EaiK) ~ EpiK)] [E,,{K)-E^{K 



X Tr 



[Pa(i^) 






dKa 



dKb 



,,,,d^(K) ,^^,d^(K) ,^^, 



dKa 



dKb 






dKa 



dKb 






dK, 
d^(K) 
dK, 

dMK) 
dK, 

d.MK) 



dKa 



dKt 



dK, 



-{b^c) , 



(40) 



The proof of (38) is given in Appendix B. From (38) we 
get convenient relations for the derivatives of Ya^ab{K), 
see Appendix B, 



d i 1 



Tr 



-Tr 



1 



d^^iK) 



1 



d^iK) 



(C - MK)) dKadKb (C - MKW dK, 
1 d^iK) 1 d^{K) 



iC~MK)) 



dKa 

1 



iC-MK)) 

d^{K)- 



(C-^(A'))2 dK, 



dKb 

(b^c) 



(39) 



see Appendix B. From both, (39) and (40), we can easily 
check the divergence condition (33). 

In the following sections wc shall use (37) and (40) 
to calculate numerically the geometric flux densities and 
their derivatives for metastable H atoms. We will be espe- 
cially interested in the flux densities in three-dimensional 
subspaces of K space. We will, for instance, consider the 
cases where the electric field £ is kept constant and only 
a magnetic field B varies or vice versa. The geometric flux 
densities (29), (37) are then equivalent to three-dimen- 
sional complex vector fields. Indeed, let us consider the 
case that only three components of K , Ka^ , A'^^ and A'aj , 
are varied. The vectors 



(41) 



span the effective parameter space which is now three di- 
mensional. In (41) we multiply the Ka^ with constants 
1/ Hi which, in the following, will be chosen conveniently. 
We shall, for instance, always choose the k^ such that the 
Li have the same dimension for i = 1, 2, 3. Wc define the 
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geometric flux-density vectors in L space as 



We denote by TZ its representation in the n = 2 subspace 
of the hydrogen atom. We have then 



Ji'HL) = 






n^{S,B)n-^ ^^iRS,RB) . 



(48) 



(42) 



where i,j, k g {1, 2, 3}. The curve C (25) and the surface 
J-" (27) live now in L space. With the ordinary surface 
element in L space 



d/f 



-e^Jk dLj A dLk 



we get for the geometric phase (28) 

7a(C)=^jW(i:)d/^ 

From (33) we find that 

div4^)(Z:)=0 



(43) 



(44) 



(45) 



AL) 



wherever (11) holds. That is, the vector fields Js, can 
have sources or sinks only at the points where the com- 
plex eigenvalues (10) of ^{K) become degenerate. More 

precisely, we see from (37) that Ja can have such sin- 
gularities only where Ea[K{L)) becomes degenerate with 
another eigenvalue Ep{K{L)) (/3 ^ a). This is, of course, 
well known [1]. From (36) we can also calculate the curl 



ofJ^) 



(rotjW(i))^ 



^ i(L 



^n'(i) 



dKa^ 



Ya,,a^a,{K{L))K,K] 



(46) 



Knowledge of both, div J4 and rot Jq , will allow us an 
easy understanding of the behaviour of the geometric flux 
density vectors for concrete cases in Section 5 below. 



4 Structure of phases and flux densities from 
rotational invariance 

In this section we shall discuss what we can learn from 
rotational invariance about the geometric phases and flux 
densities for the metastable hydrogen states. 



4.1 Proper rotations 

We consider the mass matrix ^^{£{t),B{t)) of (2). Let 
R G S0(3) be a proper rotation 



IX . Xi r J^ij Xj , 



R^{Rij), deti?=l 



(47) 



This shows that ^(5, B) and ^{R£, RB) have the same 
set of eigenvalues. Since we have assumed non-degeneracy 
of the eigenvalues of ^(£, B), see (11), the same holds for 
^^{R£,RB). Moreover, S0(3) is a continuous and con- 
nected group, therefore the numbering of the eigenvalues, 
as explained in Appendix A, cannot change with R. Thus, 
we get 



EaiR€,RB) =Eai£,B) 



For the resolvent, cf. (15) with n = 1, wc find 



(49) 



n{c-^(e,B)) ^nr^ = {c-^{R£,RB)) \ (50) 
J2{C-Ec,{£,B)y'ni>^{£,B)n-' 

^Y{C-EciR£,RB)y^¥a{R£,RB) . (51) 



With (49) we get from (51) 



nFa{£,B)Er^ = 'Pa{R£,RB) . 



(52) 



Note that for the states themselves we can only conclude 
that \a,R£,RB) and ]Z\a,£,B) must be equal up to a 
phase factor. 

With the identification of {£, B) and K of (5) we can 
now decompose the 6x6 flux density matrices Ya^ab (29) 
into 3x3 submatrices corresponding to the £ and B and 
mixed £, B differential forms (sec Appendix C of [12]). We 
have with a e /, the index of a metastable state. 



Ya,abiK) dKa A dKb = ll%(5, B) d£j A d£k 
+ I^^X(£^B)dB,AdBu 



+ r^^i^>i£,B)d£,AdBk, (53) 



where 



^L%i£^^)+^L%i£^^) = o^ 



a.kj ^ 



^n(^.»)+^n^,(^.») = o- 



a.kj *> 



(54) 
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From (2), (29), (37), (38), (53) and (54) wc obtain 
X(^)(5,B) = (45,(5,6)), 

X Tr[P„(5,B)^j.P^(£,B)^, - (j o k)] 



i 1 
227ri 



dCTr 



1 



C-^(£,s) 
1 



""-'{C-Me.B)?-' 



(55) 






I3=ta 



Tr[P„(5,B)M.PM'^,eV, - (j ^ k)] 



\^J d^Tr 



X /x. 



1 



1 



-^■(C-^(5,e)7^'-^ 

j(^'«)(£,B) = (xSf(5,B)), 

xTr[P,(5,B)^j-P^(£,B)/£^, 

-P„(5,B)/£^P^(5,B)^^-] 



(56) 



. 1 

^27ri 



dCTr 



X D, 



1 



C-^(5,B) 

1 

(C-^(£,B))2^fe 



Mt 



(57) 



using 



Wc set 



n-^D^n 



R,kD 



Jk^k ' 



^ V E: = ^jfcA^, 



RI^^\e,B)R'^ =I^^'>{Re,RB) , 
i?li^'^) (£, B) E^ ^ i^f ,e) (^^^ ^g) 



(60) 



(61) 



4.2 Improper rotations 



For improper rotations, R e 0(3) with dcti? = —1, we 
have no invariance due to the PV term E j^^'V in the 
mass matrix (3). In fact, we are especially interested in 
PV effects coming from this term. In this subsection we 
shall decompose the geometric flux densities (53) to (58) 
into PC and PV parts. 

It is clearly sufficient to consider just one improper 
rotation, the parity transformation 



P : X ^> P X = X 



(62) 



In the n ~ 2 subspace of the hydrogen atom this is rep- 
resented by a matrix P which transforms the electric and 
magnetic dipole operators as follows 






(63) 



The mass matrix ^ ^ (3) is, of course, not invariant under 
P and we have 



P^o^"'--^o, 


(64) 


PJ(pyV^^ = -^pv , 


(65) 


V^^V-^ = V{J^a + S-^Pv)^-^ 




= .£^-5.£^y . 


(66) 



As explained in general in (41) ff. we introduce the geomet- 
ric flux- dens 
components 



ric flux-density vectors J a {£,B) and Ja {£,B) with 






(58) 



see also Appendix C of [12]. 

From (53) to (57) we get the decomposition of the 6x6 
matrix {Ya,ab) in terms of 3 x 3 submatrices 



r(£) 



{Yo.,ab) 



-hd^^^Y 






T-(8) 



(59) 



The rotational properties of Iq , Iq and Xa ' are 
now easily obtained from (48) to (52) and (55) to (57) 



Clearly, since 5 w 7.57 x 10~^^ is very small (see Appendix 
A), it is useful to consider the case where it is set to zero, 
that is, where parity is conserved. We denote the quan- 
tities corresponding to this case by ^^°', Ea , Pa etc. 
We have then from (2) and (3) 

M.^°\S,B)=^o^S.-£-fJ^-B , (67) 

^{S, B) = ^("' (5, B) + (5^pv • (68) 

For the PC quantities wc find from (63), (64) and (67) 

V^^°\£,B)r-^ =^(°'(-£:,B) , (69) 

4")(£,S) = 4«)(-£,B), (70) 

'P'¥^^\£,B)V-^ =p(,°)(-£,B) . (71) 

Here (70) and (71) need some further discussion; see Ap- 
pendix C. 
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Due to time reversal (T) invariance^ and eondition (11) 
Ea{£, 13) gets no contribution linear in d, that is, linear in 
the PV term S ^py\ see [51,52]. Neglecting higher order 
terms in S wc have, therefore. 



and analogous expressions for Xa and Xa ■ These, 

the derivation of (81), and further useful expressions for 
PV fluxes, are given in Appendix C. 



£;„(£, B) -i?W (5, B) = E(^H-£,B) 



(72) 4.3 Expansions for flux densities 



Rotational invariance (49) implies then that Ea{S, B) must 
be a function of the P-even invariants one can form from 
£ and B: 



E^{£,B) = E^{£\B\{£-Bf) 



(73) 



In Subsection 4.3 below we present the analogous anal- 
ysis in terms of invariants for the geometric flux densities 
(55) to (58). For the case of no P-violation we easily find 

from (63) and (69) that zi^' {£, B) and x'^^ {£, B) must be 

even, whereas Xa ' {£,B) must be odd under £ — >• —£. 
This allows us to define generally, without any expansion 
in S, the PC and PV parts of the geometric flux densities 
as follows: 



l(f)PC(£,B) = i [xi^\£,B)+xi^\-£,B) 
I^)PV(£,B) = i [xi^\£,B)-Xi^\-£,B) 
J(«)PC(£,B) = I [if )(5,6) +If )(-5,B) 
If )P^(5,B) = I [if )(£,6) -If )(-5,B) 

If ■^)pc(£,B) = 1 [if ^)(^,e) -If -^H-^,©)' 



If ■«)PV(£,B) = ^ [if «)(£,B)+If ^^H-'^,^) 



These combine to 

I„(£,B)=lPC(£,B) 



ir(£,B) 



(74) 
(75) 
(76) 
(77) 

(78) 
(79) 

(80) 



for Iq = Iq , Iq and if. 

Using now the expansion in the small PV parameter 
S up to linear order we get for the PC fluxes exactly the 
expressions (55) to (57) but with ^{£,B), ]Pa,i3{£,B), 
Ea{£, B) replaced by the corresponding quantities for 5 = 
0, that is, ;y£}-^'{£,B) etc. For the PV fluxes we have to 
expand the expressions (55) to (57) up to linear order in 
5. This is easily done and leads to 



"a.Jk 






-£>, 



dCTr 



C-^^'"i£,B) 



-PV 



£,B)-'{C-^^'>\£,B)) 



iDk 



U ^ k) 

(81) 



^ Here we disregard the T-violating complex phase in the 
Cabibbo-Kobayashi-Maskawa Matrix. This is justified since we 
are dealing with a flavour diagonal process. 



We can now write down the expansions for the geometric 
flux densities (55) to (57) following from rotational invari- 
ance and the parity transformation properties. In these 
expansions we encounter invariant functions 

g'^ = g^{£\B\{£-B)'), 

h^ = h^{£\B^,{£-B)^) , 

r = l,...,15, (82) 

which are, in general, complex valued. Our notation is 
such that the terms with the g" are the parity conserving 
(PC) ones, the terms with the /i" the parity violating (PV) 
ones. We find the following: 



r(£) 



1 



t(£), 



X^>,i£,B)^-e,,iJ^>i£,B), 

jf)(£,B) = jf)PC(£,B) + jf)PV(5,B) 

= £[(£■ B) g? + h<(\+B [gq + {£-B) h^] 
+ {£xB)[{£-B)gS + h^], (83) 



^f(^.e)=i: 



£,B)PC 

a.jk 



{£,B)+X\ 



(£,B)PV 
ajk 



i£,B) 



(84) 



= l^jki {Si [g^ + {£-B) h^] + Bi [{£ ■ B) gt + ht\ 
^{£xB\{gt^{£-B)l4\]^ 
+ 5,fc[(5-B),9? + /i?] 

^{£,£u-\&,k£''){{£-B)gt-rl4\ 
+ {B,Bu-\&,kB'){{£-B)gt^h%\ 

+ {£,Bk + £uB, - \b,u {£ ■ B)) [.g?o + (^ • ») K^\ 
+ {£, {£ X B)k + £k {£ X B)j\ [{£ ■ B) g'^, + h^,] 
+ [B, {£ X B)k + Sfe (£ X B),\ [ff?2 + (^ • 6) h^^ 

^SU'^,S) = i6,fc,jl^)(5,S), 

J'^^'^ {£, B) = Ji^^''''{£, B) + Jl^^''^{£, B) 

^£[{£-B)g^, + K^] 

+ B[g^^ + {£-B)K,] 

+ {£xB)[{£-B)g^^ + h%]. (85) 

In the next section we shall discuss specific examples 
where the general expansions (83) to (85) will prove to be 
very useful. 

5 Specific Cases 

In this section we shall illustrate the structures of the ge- 
ometric flux densities for specific cases. First, we analyti- 
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cally derive the geometric flux densities Ja {£■ = 0, 3) in 
magnetic field space for vanishing electric field and com- 
pare the results with numerical calculations. In that case 
the flux densities give rise only to P-conserving geomet- 
ric phases. We then analyse the structure of Ja {£,3) 
in electric field space together with a constant magnetic 
field B = S3 63. Here, we will obtain P-violating geometric 
phases. As a further example we investigate the geometric 
flux densities in the mixed parameter space of £1 , £3 and 
B3 together with a constant magnetic field 6 = ^262. 



5.1 A magnetic field B and £ = 



(B), 



Starting from the general expression (85) for Ja (5,B), 
we immediately find for Ja {£ = 0, 3) 

j(«)(0,B)=Bg?4 (86) 

where ^"4 = 3^4(6 ) only depends on the modulus squared 
of the magnetic field. 

Because of (45) we have 



Vb-J(^)(O,B)=0 

for 6 = B V 0. Inserting (86) in (87) we get 

= Vb • {3gUb)) = igUb) + 2bdbgUb) 
with the solution 

<4(6)=a"fe-i=a"|B|-3 



(87) 



(89) 



where a" is an integration constant. Therefore, from the 
rotational invariance arguments above we find already 



j(^)(0,B) = a° 



3 

IBP 



(90) 



This is the field of a Dirac monopole of strength a" at 
3 = 0. A detailed calculation of a" for the 2S states is 
presented in Appendix D. The resulting flux-density vec- 
tor field is real and P-conserving. It reads 






for Q = 9 



JL^Ho,3)={ -^ , fora^ll, 
, for a = 10, 12 



(91) 



We recall that the states with labels a = 9, 10, 11, and 12 
are connected to the 2S states with {F, F3) = (1, 1), (1, 0), 
(1,-1), and (0,0), respectively; see Appendix A, Table 2 
ff. 

Comparing this exact result with numerical calcula- 
tions, we can extract an estimate for the error of Ja in 
parameter spaces other than that of the magnetic field. For 
lO'' equidistant grid points in a cubic parameter space vol- 
ume [—1 mT, 1 mT]^, at which the vectors j}]^ (0, 3) are 



0.5 



o 
2Q 



K) 



-0.5 



_ • • ' / \ ^ ■ • _ 



-0.5 





Bi/Bo 



0.5 



Fig. 3. Visualisation of the P-conserving flux-density vector 



field ji^'^^(0,B) (94) for a = 11 and 77 
field parameter space at B3 =0. 



10 in magnetic 



evaluated, we obtain numerically the following deviations 
from the vector field structure given in (91): 



j{f (0,B)|-|Bp-l| <5x 10-12 , 
j[f\o,3) 3 



Jif\0,3) 



\3\ 



< 8 X 10 



-13 



(92) 



Since the data type used for the numerical calculations has 
a precision of approximately 16 digits, we find the results 
(92) to be in good agreement with the analytic expression 
(91). 

Figure 3 illustrates the numerical results for Ja (0, 3) 
with a = 11 in the B3 ~ plane. 

Here and in the following we find it convenient to plot 
dimensionless quantities. Therefore, we choose reference 
values for the electric and magnetic field strengths 



£0 = 1 V/cm 
60 = 1 mT . 



(93) 



We label the axes in parameter space by Bi/Bo and plot 
the vectors 



Jl^) 



V 



J^^^Bl. 



(94) 



Here 77 is a rescaling parameter chosen such as to bring 
the vectors in the plots to a convenient length scale. The 
dimensionless geometric phases, see (26), (28), and (44), 
are then given by the flux of this dimensionless vector field 
(94) through a surface in this space of Bi/Bo and divided 
by 77. 
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For any curve C ~ dJ- in B space we get for the geo- 
metric phases from (44) and (91) 



7.(C)=/ Jf)(B)d/(«) 



fie, 


for a = 9 , 


^c, 


for a = 11 , 


. 


for a = 10, 12 



(95) 



Here f2c is the sohd angle spanned by the curve C. This 
is in accord with the expectation for a spin 1 system; see 
[2]. 



5.2 An electric field £ together with a constant 
magnetic field 

We now consider the case of geometric flux densities in 
electric field space with a constant magnetic field B = 
Bs 63 with S3 > 0. Here we find from (83) 



7(^)^^(5,^363) =5^:3^3 3?+ S3 



f 2^3^3.9? 

-£i£3B3g^ 
9^ 



J(^)PV(£, S3 63) = £K + Bs 63 £3^3 h^ 

+ (£•261 -£163)^3/1^ 

/ £,h- 
= £hf + B3\ -£i h^ 



(96) 



(97) 



where gf2 3 ^^'^ '*?2 3 ™^3,y in general depend on 5 , S| 
and (£3^3)^ and a G {9, 10, 11, 12}. In our case B3 is con- 
stant. We have, therefore. 



a a / c»2 c2\ 

9i =9i (^T^Ss) 
hf^hf{£'r,£i] 
(z = 1,2,3) 



with 



C- f-p - — Oi 



fl- 



(98) 



(99) 



In the following we give the results of numerical eval- 
uations of the PC and PV flux-density vectors (96) and 
(97), respectively. We split the PV vectors into the con- 
tributions from the nuclear-spin independent (i = 1) and 
dependent (i = 2) PV interactions 



j(£)PV ^ j(£)PVi _^_ j(£)PV2 



(100) 



Here the J(^)p^» are defined as in (81) but with S replaced 

by Si and ^py replaced by ^py {i = 1,2); see (A. 5) and 
(A. 6). Again we shall plot dimensionless vectors 



j(^)P^(5,S3e3) 
j(^)P^-(5,S3e3) 



V'a.i 



j(^)PC(5,S3e3)£o', (101) 
j(^)P^-(5,S3e3)£o' (102) 



0.5 






-0.5 






-0.5 



0.5 



^2/^0 



Fig. 4. Visualisation of the real part of the P-conserving flux- 
density vector field Jq ' {£,B-j, e-j) (101) in electric field pa- 
rameter space at £1 = 1 V/cm. A constant magnetic field with 
Bs = ImT is applied. The scaling factor in (101) is chosen as 
r?9 = 2.5 X 10'*. 



with £{) from (93) and rja and r]'^ ^ conveniently chosen. 

As an example of a PC fiux-density vector field we 
present in Figure 4 the results of a numerical calculation 

(£)PC, 



of Jq'^'''^{£, B3 63) for B3 = 1 mT. We recall that the state 
with a = 9 is connected to the 2S state with {F, F3) = 
(1, 1); see Appendix A. We plot the dimensionless vectors 
(101) with the scaling factor chosen as 779 = 2.5 x 10^. 
Comparing with (96) we see that here the dominant term 
is the one proportional to g2'- 



(103) 



jf)P^(5, 6363) « S3 g2'('^T, ^1)63 

with 32(^X1 '^3) being practically constant. 

The numerical results shown in Figure 4 reveal a large 
sensitivity of the flux-density vector field Jg to the 

normalisation of the dipole operator D. Indeed, suppose 
that we make in our calculations the replacement 



D^XD 



(104) 



with A a positive real constant. From the mass matrix 
(2) and from Zg ' in (55) we find the following scaling 



behaviour for J^ 



jf)^^(£,S3e3) 



r(£)PC 



\D 



A^ jr^""(A£,63e3) 



(105) 



Since our calculations show that J, 



:£)PC 



is practically con- 



stant for the range of fields explored here we get 



J^ (5,i33e3) « A J^ ' (5,^363) 



(106) 
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Fig. 5. Visualisation of the real part of the P-conserving flux- Fig. 6. Visualisation of the real part of the nuclear spin depen- 



density vector field Jg ' (£,83 63) (101) in electric field pa- 
rameter space at £^1 = 1 V/cm. A constant magnetic field with 
B3 — 1 fiT is applied. The scaling factor in (101) is chosen as 
rj9 = 450. 



dent P-violating flux-density vector fleld J, 



{£)PV2 



{£,B) (102) 



in electric field parameter space at £3 = 1 V/cm with an ad- 
ditional constant magnetic field B — Bse^, Bs — IfiT. The 
scaling factor in (102) is chosen as 779 2 = 400/(52. 



Therefore, measurements of Jg {£, B3 63) for the setup 
considered here arc highly sensitive to deviations of the 
normalisation of D^ from the standard expression as given 
in Table 4. 

As an example we calculate the PC geometric phase 
for the following path in £ space 



1.0 + 0.5 X sin(200 7rz) 

0.5 X cos(200 7rz) 

1.0 



That is, we consider a path circling 100 times in the 
Si — £2 plane which is orthogonal to the flux direction for 
a = 9; see (103) and Figure 4. For the constant magnetic 
field chosen there, 83 = ! mT, we obtain for a = 9 





79(C) = 5.89 X 10""^ - 5.25 x 10"^ i 



Decreasing the constant magnetic field B3 we get larger 
geometric phases. Numerical results of the flux-density 
vector field Jg ' {£,83 e^) for B3 = 1 ^T are presented 
in Figure 5. There, the scaling factor is chosen as ijg = 450. 
For the curve (107), S3 = 1 fiT, and a = 9 we obtain 



79(C) =0.0340-0.005961 



(109) 



As an example of a PV flux-density vector field wc 
present the numerical results for Jg ^ in Figure 6. We 
find both, the real and the imaginary part of 
Jg ^(5,^363) to represent a flow, circulating around 
the e3-axis, with vanishing third component. That is, in 



(97) for a = 9 the terms involving h\ and /12 , respectively, 
come out numerically to be negligible compared to the 
terms involving /13. We may hence write 



r(£)PV2 



£2 

{£,B3e3)~B3hl{£l,£l) \ -fi 





(110) 



This is corroborated by numerical studies. For \£j \ < £q 
(j = 1, 2, 3) with £q from (93) wc find 



|Ree3jf)^^ 



I Re J 



(£)PV2| 
9 I 



< 6 X 10" 



and 



(108) |gi Re ei jf ^^""^ + £2 Re e^ jf ^^""^ | 



I Re J 



(£)PV2| 



VW+^: 



(111) 

< 2.5 X 10~" . 

(112) 



Results similar to (111) and (112) hold for the imagi- 
nary part fmJg ' ^(£,^3 63) and for Jg ' ^(£,^363). 

Thus, also Jg ^ has to a good approximation the struc- 
ture (110). 

The antisymmetry of Re Jg ^ (5 , iJs 63 ) with respect 
to 5 — > —£, see (97), is confirmed numerically at the same 
level of accuracy. We find 



|Rejf^P''^(£,63e3)| 



|Rcjf)P^^(-£,S3e3 



< 2.5 X 10" 



(113) 
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and 



|Re Jf ^^"'^ (5, S3 63) + Re jf ^^""^ {-£, S3 63) 



|Rcjf)P^^(5,63e3)| 



< 2.5 X 10" 



(114) 



Similar numerical results are also obtained for 



Im J( 



(£)PV2 



9 



{£, B3 63) and for the real and imaginary parts 



ofjf)P^^(£,63e3) 



5.3 The mixed parameter space of fi, £3, B3 together 
with a constant magnetic field 

In our last example we discuss the parameter space spanned 
by the vectors 



(115) 



see (41) ff. and (93). In addition we assume the presence 
of a constant magnetic field S = 1 jiT 62 . The values cho- 
sen for £q and Bq in (93) and (115) should represent the 
typical range of electric and magnetic field variations, re- 
spectively, for a given experiment. Our choice here is mo- 
tivated by the discussion of the longitudinal spin echo ex- 
periments in [14]. For other experiments different choices 
of £0 and Bo will be appropriate. 

From (41), (42), and (55) to (59) we find 





J(^)(5,B) = 



£oBo 











£oBo 














~li'if{£,B) 
2li%{£,B) 



(116) 



Here again a G {9, 10, 11, 12} labels the metastable state, 
see Table 2 of Appendix A, to which the geometric flux 
density corresponds. Inserting in (116) the expressions 
from (83) and (84) we get for the PC and PV parts of 

ji^^(5,B) 




B3£3g?+B3£i9^ 

B3£3gS+B3£i92 
~gi+£i£3~9i 



hi + £i£:3 h^ 

h^ + £i£3 K 

B3£3 k + ^83^1 K . 



(117) 
(118) 



where the functions gf and hf, i £ {1,...,6}, may in 
general depend on £f, f| and B^; see Appendix D. 

For ease of graphical presentation we shall multiply 
the Ja with scaling factors rj. Thus we define 

J^^)P^(£,B) = ^„J(^)P^(5,B), (119) 

(i = l,2). (120) 
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Fig. 7. Visualisation of the 1 and 2 components of the real part 
of the P-conserving flux-density vector field Jg (f,B)(119) 

at Z?3 — 143 /iT. The scaling factor is chosen as r/g — 2.5 x 10*. 
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Fig. 8. Visualisation of the 1 and 2 components of the 
imaginary part of the P-conserving flux-density vector field 
j^i,)pc^g g^ (119) at B-i = 143 AtT. The scaling factor is cho- 
sen as 779 = 2.7 X 10^. 



In the Figures 7 to 11 we illustrate the results of our 
numerical calculations of (119) and (120) for the geomet- 
ric flux densities of the state with a = 9, that is, the 
state connected with the 2S, (F, F3) = (1,1) state; see 
Appendix A, Table 2. For 30'^ grid points in the parame- 



ter space volume \—£o^ ^0] 



-0.2^0, 0.2i3o] we find that 
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Fig. 9. Visualisation of the 1 and 2 components of the real part 
of the nuclear spin dependent P-violating flux-density vector 



{-L)PVi 



Fig. 11. Visualisation of the 1 and 3 components of the real 
part of the nuclear spin dependent P-violating flux-density 



field J, 

is chosen as 779 2 = 10''/J2- 



(e,B) (120) at S3 = 143 ^T. The scaling factor vector field J, 



(i)PV2 
9 



(£-,16) (120) at £z = 0. In order to re- 
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solve the structure more clearly the scaling factor is chosen as 
r,^,2 = 2xl05|S3/So|V^2- 



in the E\ — £3 plane as representative vector field struc- 
tures. Investigating the dependencies of Jg ' (£, IS) on £1, 
£3 and S3 we find our numerical results to be consistent 
with the analytic structures in (117) and (118). We recall 



from (45) that J< 



(i)PC 



is divergence free and, thus, gen- 



crated by a non- vanishing vortex distribution rot J^ 
We have calculated this quantity numerically using (46) 
to be in agreement with a direct eval- 



and find rot Ji 



(Z,)PC 



9 



(i)PC 



-0.5 





^1/^0 



0.5 



nation using a fit function representing J, 

Examples of PV flux-density vector fields are shown in 
Figures 9 to 11. Again, the full three-dimensional vector 
fields are divergence free and thus represent flows without 
sources or sinks. Note that in Figure 11 we have chosen 
a scaling factor depending on S3 since the magnitude of 



the vector J, 
S3 = 0. 



i)PV2 



(£, S) shows a strong increase towards 



Fig. 10. Visualisation of the 1 and 2 components of the real 
part of the nuclear spin dependent P-violating flux-density vec- 
tor field J^^^^^2 (f , JB) (120) at Bz = 86 ^T. The scaling factor 
is chosen as 779^2 = 4 x 10^/(52. 



numerically 

|Rce3jf'P^|<0.2[(Reeijf)P^)2 

+ (Ree,jf)P^)2]\ (121) 

An analogous relation holds for Im Jg (6,13). This jus- 
tifles the presentation of these flux-density vector fields 



The emergence of non- vanishing imaginary parts of the 
flux-density vector fields, for example, in the parameter 
space considered in the present section, see Figure 8, leads 
to an interesting phenomenon. Let us consider a closed 
curve C, being the boundary of a surface J-', in L space 
(115). We parametrise C as 



C: z^ L{z) , 
0< z<l , 

Z(0)=i(l), 



(122) 



cf. (25), (27). We suppose that over a time interval [0,T] 
this curve in parameter space is run through by the system 
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in the following way: 



0<t<T , 

C:t^ L{z{t)) . 



(123) 



We consider the adiabatic limit where T becomes very 
large. We shall also consider that the system is run through 
the curve C in the reverse direction; 



t -^ z{t) = 

0<t<T, 

C: t^ L{z{t)) 



T-t 



T 



1 - z{t) , 



(124) 



Suppose now that we have at t = the atom in the 
initial state i/'a(0), see (21). We change the parameters L 
along the curve C as in (123). From (21) to (23) we find 
the decrease of the norm of the state at time T to be 

1^! (0) I ^ = cxp [+2 Im y^ (T) - 2 Im 7„ (T)] . (125) 

Here 2 Im (pa (T) and 2 Im 7^ (T) are the contributions due 
to the dynamic and geometric phases, respectively, 

2l-mipa{T)=T2lm dz Ea{L{z)) , 
Jo 

2Im7,(r) = 2Im /" jl^\L)df^ ; (126) 

see (26) and (44). From (125) we can define an effective 
decay rate for the state a under the above conditions as 

r„,eff(C,r) = ^ [-2Im(^„(T) + 2Im7„(r)] 

= -2Im / AzE„{L{z)) 
Jq 

+ |;Im^jW(i)d/^- (127) 

Note that this effective decay rate depends, of course, on 
the curve C and that the geometric contribution is sup- 
pressed by a factor 1/T relative to the dynamic contribu- 
tion. From (125) and (127) we get for the decrease of the 
norm of the state 



^-^^=exphr„,eff(C,r)T] 



(128) 



Now we start again with the state ipa{0) at time t — 
but we change the parameters L along the reverse curve 
C (124). It is easy to see that the dynamic term in (127) 
does not change whereas the geometric term changes sign. 



r„,eff(C,r) = -2Im / dzE,,{L{z)) 



T 



Im / Ji'^>{L)df 



(129) 





B3 [/iT] A 
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/£^ [V 



£^ [V/cm] 

Fig. 12. The curve C (130) in the parameter space £1, £^, B3. 
The circle is run through twice. 



Thus, the effective decay rate depends on the geometry 
and reversing the sense of the running through our closed 
curve in parameter space changes the sign of the geometric 
part. 

As a concrete example we choose a constant magnetic 
field B2 62 with B2 = I fJ-T and the following curve in L 
space 



f£iiz)/£o 

C: z^L{z)^ £z{z)/£o 
\B3iz)/Bo^ 

£i{z) =0.8 V/cm , 

£3(2) = 0.5 X cos(4 7rz) V/cm , 

iJg (z) = [2 + 2 X sin(4 TT z)] ^T , 

0< z< 1 ; 



(130) 



see Figure 12. We suppose as in (123) that C is run through 
in a time T with z(t) = t/T. In this time the path in pa- 
rameter space makes, according to (130), two loops. With 
T = 1 ms we can meet the adiabaticity requirements as 
spelt out in [14] and in Equations (30), (31), and (37) of 
[12]. The essential requirement here is that the frequency 
1/ = 2/T of the external field variation in (130) must be 
much less than the transition frequencies AE/h between 
the Zeeman levels for a = 9, 10, 11. For the external field 
of order 1 fiT we get AE/h ~ 10 kHz which gives the re- 
quirement 



1.== - < lOkHz , 
r> 0.2 ms . 



(131) 



J" 



Calculating now the contributions to the effective decay 
rate (127) we find for the state a ~ 9 which is connected 
to the 2S state with (-F, F3) = (1, 1), see Appendix A, the 
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following 



2Ini / dzE^mz)) = 1935.2 s"^ , (132) 

;Im / jf^(i)d/^ = -1.8s-i . (133) 



This leads to 

r9,eff(C,r) = 1933.4 s-i . (134) 

For the reverse curve C wc find, instead, 

Tg^off (C,T) = 1937.0 s-i . (135) 

Thus, under the above conditions the effective decay rates 
(134) and (135) differ by 1.9 %o and the corresponding 
decreases of the norms (128) by 3.6 %o. We emphasise that 
this difference has its origin in the geometric phase. 

To give an example of a PV geometric phase wc con- 
sider the following curve 

( £i{z)/£q 

C: z^L{z)^ £z{z)/8a 
\Bsiz)/Bo^ 

£i{z) = 0, 

Ssiz) — So sin(2 7rz) , 

Bsiz) = 0.16o cos(2 7rz) , 

< z < 1 . (136) 

For this curve the PC geometric phases vanish due to the 
antisymmetry of ei ■ Ja (5, S3 63) under {£1,8:^) — ^ 
(— fi, —£3); see (117). Thus, we get here 

Numerically wc find for a — 9 and S2 = 1 mT 

7^^i (C) = (0.00467 - 0.000457 i) Si , (138) 

7^^" (C) = (0.0942 - 0.00421 i) S2 . (139) 

Assuming now that we can circle the curve C N times 
we get as geometric phase N'-fg{C'). Thus, the number of 
circlings acts as an enhancement factor for the small weak 
interaction effects in hydrogen. For N = 10^, for example, 
we obtain 

7V7^^(C') = (46.7 - 4.57 i) ^i + (942 - 42.1 i) 62 ■ (140) 

With Si 2 from Table 1 this gives a phase of the order of 
10-9. 



6 Conclusions and outlook 

In this article wc have discussed the geometric phases and 
flux densities of the metastable states of hydrogen with 
principal quantum number n = 2 in the presence of ex- 
ternal electric and magnetic fields. We have provided ex- 
pressions for the flux densities and their derivatives suited 



for an investigation of their general structure. This was 
achieved with the help of representations of the flux den- 
sities as complex integrals. For these integrals extensive 
use was made of resolvent methods, and the results turned 
out to be quite simple and easy to handle. Furthermore, 
employing proper and improper rotations we derived the 
general structure of the flux densities for the metastable 
states. We also obtained expansions of the flux densities 
in terms of P-conserving and P-violating contributions. 
The flux densities can be visualised in the case of three- 
dimensional parameter spaces as vector flclds. We gave 
three examples of parameter spaces for which we com- 
pared analytical with numerical calculations. The results 
are consistent regarding the employed numerical precision. 
For vanishing electric field the flux densities in magnetic 
field space are real and P-conserving - and so are the cor- 
responding geometric phases. In this case the flux-density 
vector field can be computed analytically and turns out 
to be the field of a Dirac monopole sitting a.t B — 0. 
In parameter spaces of both electric and magnetic field 
components the flux densities exhibit a rich structure of 
P-conserving as well as P-violating vector fields including 
both real and imaginary parts. In the general structures 
of the fiux densities of those cases we encounter functions 
which are rotationally invariant; see (83) to (85). These 
functions contain all the information on the mass matrix 
and the electric and magnetic dipole matrices which one 
can obtain from the measurement of the geometric phases 
for the system considered. 

In Section 5 we have calculated geometric phases for 
various situations and, surely, the question arises about 
their possible measurements. Present experiments can 
reach a precision in phase measurements of about 10"^ rad 
[53]. Thus, the PC phases (95), (108), (109) and the change 
of decay rate (134) and (135) should be within reach of 
these experiments. The PV phases for hydrogen certainly 
need further theoretical and experimental efforts to bring 
them to a practically measurable level. 

The general representations for the fiux densities as 
complex integrals are, of course, easily transferred to other 
atomic systems having stable or metastable states, for in- 
stance, to the states of n = 2 helium. Similarly, the anal- 
ysis of the consequences of rotational invariance and of P 
violation in Section 4 goes through unchanged for other 
atomic systems. The only requirement here is that the 
coupling of the system to the external electric and mag- 
netic fields can be described as in (2) with electric and 
magnetic dipole matrices D^ and /x, respectively. We note 
that measuring geometric phases for such systems can give 
valuable information on their atomic matrix elements. We 
have found for hydrogen that, on the one hand, there is the 
case of a pure magnetic field where the flux densities and 
geometric phases are simply given by the geometry of the 
path run through in parameter space; see Section 5.1. On 
the other hand, a large sensitivity on the electric dipole 
matrix element was found for the flux densities in elec- 
tric field space with the presence of a constant magnetic 
field in Section 5.2. We also discussed the change of the 
effective decay rates for geometric reasons in Section 5.3. 
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All these phenomena should also occur for other atomic 
systems. Measurements of geometric phases could be an 
important testing ground for theoretical calculations of 
wave functions and matrix elements for atomic systems, 
for instance, for He. 
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Appendices 

A The n = 2 states of hydrogen 

In this appendix we collect the numerical values for the 
quantities entering our calculations for the hydrogen states 
with principal quantum number n = 2. We specify our 
numbering scheme for these states. The expressions for 
the mass matrix at zero external fields and for the electric 
and the magnetic dipole operators are given in Appendix 
E. 

In Table 1 we present for }H, where the nuclear spin 
is / = 1/2, the numerical values for the weak charges 

Q^ , >f = 1,2, the quantities Aq + Aq, the Lamb shift 
L = S(2Si/2) -^(2Pi/2), the fine structure splitting A = 
S(2P3/2) — £^(2Pi/2), and the ground state hyperfine split- 
ting energy A. We have A = £'(lSi/2, F = 1) - -B(lSi/2, 
F = 0) for hydrogen. We define the weak charges as in 
Section 2 of [12] which gives for the proton in the stan- 
dard model (SM): 






= 1-4 sin^ dw , 
X {Au + Au- Ad- 



Ad - As - As) 



(A.1) 



Here i^vv is the weak mixing angle and Aq + Aq denotes 
the total polarisation of the proton carried by the quarks 
and antiquarks of species q {q ~ u,d,s). Note that in 
[12] and [47] we adhered to the then usual notation of 
Aq for what is now denoted as Aq + Aq. The quantity 
Au + Au — Ad — Ad is related to the ratio qa/ gv from 
neutron /3 decay: 



Au — Ad — Ad = —gA/fJv 



(A.2) 



The numerical value given in Table 1 is from [27]. The total 
polarisation of the proton carried by strange quarks. As + 
Zis, is still only poorly known experimentally. One finds 
values of —0.12 to very small and positive ones quoted 



in recent papers; see for instance [54-58]. Therefore, we 
assume for our purposes 



-0.12 < As-|-Z\s< 



l(2) 



(A.3) 



Of course, the dependence of Q\iJ on As-\-As is, in princi- 
ple, very interesting, since this quantity can be determined 
in atomic P violation experiments with hydrogen. 

We define (see (19) of [12]) the dimensionless constants 



<5. = - 



^/iG 



q1"^ 



'w 



64:7T\/2 r'^rrie 



(><=1,2) 



(A.4) 



with Fermi's constant G, the Bohr radius rs and the elec- 
tron mass TOe . We see from Table 1 that varying As + As 
in the range (A.3) corresponds to a 10% shift in 62- Thus, 
a percent-level measurement of 62 would be most welcome 
for a clarification of the role of strange quarks for the nu- 
cleon spin. 

The mass matrix for zero external fields is given by 



+ Si^^^l + S: 



.(2) . 
-PV ' 



(A.5) 



see (22) to (24) of [12] and Table 3 in Appendix E below. 

The terms ^py and ^py correspond to the nuclear-spin 
independent and dependent PV interaction, respectively. 
As in (C.31) ff. of [12] we define 



5 = iS! 



sl)^ 



^PV == Yl J^PV 



(A.6) 



The n = 2 states of hydrogen in the absence of P 
violation and for zero external fields are denoted by |2i,/. 





IH 


Ref. 


L/h 
A/h 
A/h 


1057.8440(24) MHz 

10969.0416(48) MHz 

1420.405751768(1) MHz 


[59] 
[59] 
[60] 


5i 


0.04532(64) 
-2.78(4) X 10"" 


(11) of [12] 
(20) of [12] 


Au + Au 
-Ad- Ad 


1.2694(28) 


[27] 



As + As 

52 



-0.12 

-0.1259(18) 

7.74(11) X 10'- 



0.00 

-0.1151(16) 

7.07(10) X 10~" 



(A.3) 
(A.1)-(A.4) 



Table 1. Values of parameters for JH for numerical cal- 
culations. The weak mixing angle in the low energy limit, 
sin^ -dw = 0.23867(16), is taken from [61]. The uncertainty 

in 5 1 is dominated by the uncertainty of sin'^i^vi'- The uncer- 

(2) 
tainties quoted for Q\y and 82 for hydrogen are resulting from 

the error of the weak mixing angle. The variation of Q\^ and 
S2 with (As + As) varying in the range (A.3) is given explicitly. 
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F,F^), where L, J, F and F^ are the quantum numbers of 
the electron's orbital angular momentum, its total angular 
momentum, the total atomic angular momentum and its 
third component, respectively. 

To give the matrices ./# q, ^py , ^py , D_, and /i. ex- 
plicitly we use the following procedure. The hermitian part 



of. 



^0' 



that is, 



+ 



2 is given by the known en- 



ergy levels of the n = 2 hydrogen states; see Table 1. The 
decay matrix _r = i(^n — ^l.) needs a little discussion. 
We have (for atoms at rest), see for instance (3.13) of [47], 

{2L'j,,F',Fm'^Lj,F,F3)^27rY,{X\T\2L'j,,F',F^)* 



X 



xd{Ex-E2){X\T\2Lj,F,F3 



(A.7) 



Here, \X) denotes the decay states, the atom in a n = 1 
state plus photons, and T is the transition matrix. Ro- 
tational invariance tells us immediately that the matrix 
r^ must be diagonal in (i^, F3). We shall neglect P viola- 
tion in the decay. Then the non-diagonal matrix elements 
between S states (L = 0) and P states (L = 1) must be 
zero. The only non-diagonal matrix elements (A.7) which 
could be non-zero arc, therefore, those for L = i' = 1, 
F^l, and (J', J) = (1/2, 3/2) or (J', J) == (3/2, 1/2), re- 
spectively. But calculating these matrix elements inserting 
the usual formulae for El transitions on the r.h.s. of (A.7) 
we get zero. Thus, neglecting higher order corrections, the 
matrix r_ (A.7) is diagonal. The numerical values for the 
diagonal elements of F are taken from [37,38]. 

To calculate the matrices ^py, ^py, D, and /x we 
use the standard Coulomb wave functions for hydrogen. 
As in [47] (see Appendix B there) we use for these states 
the phase conventions of [62] except for an overall sign 
change in all radial wave functions. The matrices ^ q. 

D^ and u, in this basis are collected in Tables 3 to 5 of 

— ^ 

Appendix E (online material only). 

Now we discuss the properties, the ordering and the 
numbering of the eigcnstates of ^(5, 3) as given in (2). 

We are interested only in moderate magnetic fields, 
that is, we want to stay below the first level crossing in 
the Breit-Rabi diagram, which implies 



\B\ <53.8mT 



(A. 



Note that these crossings are only in the real part of the 
cigencncrgics; see (10). In the region (A. 8) degeneracies of 
the complex energies (10) occur for B = at arbitrary £. 
This is a consequence of time-reversal (T) invariance. For 
B = and £ 7^ we can choose the vector 63 = ^/|^| 
as quantisation axis of angular momentum. Then F^ is a 
good quantum number and time reversal invariance im- 
plies that there arc corresponding eigcnstates of ^(£, 0) 
with quantum numbers ^3 and —F^ and having the same 
complex eigenenergies. See Sections 3.3 and 3.4 of [47] for a 
proof of this result using resolvent methods. Thus, we have 
degeneracies of the complex eigenenergies of ^{E, B) (2) 
in the parameter subspace 



{(£, B) ; e arbitrary, B = 0)} 



(A.9) 



By numerical methods we checked that, at least for mod- 
erate B fields (A. 8), there are no further degeneracy points 
or regions. 



hydrogen 



\2Lj,F,Fs,£,B) 



1 


1 2^3/2,2, 


2,f,B) 


2 


|2P3/2,2, 


1,£,B) 


3 


|2P3/2,2, 


0,8, B) 


4 


|2P3/2,2, 


-1,£,B) 


5 


I2-P3/2, 2, 


-2,£,B) 


6 


1 2^3/2,1, 


1,£,B) 


7 


|2P3/2,1, 


0,£,B) 


8 


[2^3/2,1, 


-1,£,B) 


9 


|2Si/2,l, 


1,£,B) 


10 


|25i/2,l. 


0,£,B) 


11 


|2Si/2,l, 


-1,£,B) 


12 


|2Si/2,0, 


0,£,B) 


13 


!2A/2,i, 


1,£,B) 


14 


!2Pl/2,l, 


0,£,B) 


15 


l2Pl/2,l, 


-1,£,B) 


16 


|2Pi/2,0, 


0,£,B) 



Table 2. The numbering scheme for the atomic n 
of hydrogen. 



2 states 



The eigcnstates of ^ (2) for electric field £ and mag- 
netic field B equal to zero are the free 2S and 2P states. 
We write L, P, S since these states include the parity mix- 
ing due to HpY, see (1) of [12]. Thus, the eigcnstates of 
the mass matrix (2), including the PV part but with elec- 
tric and magnetic fields equal to zero, will be denoted by 
\2Lj,F,F3,£ = 0,B = 0). The corresponding states for 
the mass matrix without the PV term, that is, with ^g 
replaced by ^ p in (2), will be denoted by |2Lj, F, F3, £ = 
0, B = 0). But it is not convenient to start a numbering 
scheme at the degeneracy point {£, B) = (0, 0). Therefore, 
we consider first atoms in a constant B-field pointing in 
positive 3-direction, 



B^Be3 , B>0 



(A.IO) 



The corresponding eigcnstates, denoted by |2ijj, i^, i^3, 0, 
Be^), and the corresponding quasi projectors (12) of ^ 
in (2) are obtained from those at ;B = by continuously 
turning on B in the form (A.IO). Of course, for B > 0, 
F^ still is a good quantum number but this is no longer 
true for F. The latter is merely a label for the states. 
We now choose a reference field Bref = -BiofCs, B^^t > 0, 
below the first crossings in the Breit-Rabi diagram, for 
instance Bici = 0.05 mT. We are then at a no-degeneracy 
point and number the n = 2 states and quasi projectors 
(12) with a = 1, . . . , 16 as shown in Table 2 setting there 
{£,B) = (0, Brcf)- In the next step we consider external 
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fields of the form 



£' 





6' > , (A.ll) 



and a eontinuous path to these fields from the reference 
point (£,B) = (0,B,.cf): 



£'i\) = \£' , 

B'(A) = Sref + X{B' - Bref) , 

AG [0,1] . 



(A.12) 



Since we encounter no degeneracies for A € [0,1] the en- 
ergy eigenvalues as well as the quasi projectors are con- 
tinuous functions of A there. This allows us to carry over 
the numbering of the quasi projectors from (0, Bicf ) to all 
fields of the form (A.ll). 

Finally, we consider arbitrary fields {£, B) with B ^ 0. 
We can always find a proper rotation R such that 



R£ = £' 



RB = 



(A.13) 



with {£' ,B') of the form (A.ll). From the considerations 
of the resolvent in Section 4.1 we conclude that the eigen- 
values of ^{£,B) and ^[£' ,B') are equal. There are 
also no degeneracies here and, therefore, we can unam- 
biguously carry over the numbering of eigenvalues and 
quasi projectors from the case {£' , B') to the case (5, B). 
The labels a = 1, . . . , 16 in Table 2 for arbitrary (£, B) 
with B y^ correspond to this identification procedure 
of eigenenergies and quasi projectors. The corresponding 
eigenstates \a, £, B) oi ^M_(£^ B) arc defined as the eigen- 
states of the quasi projectors 



¥^{£,B)\a,£,B)^\a,£,B) 



(A.14) 



where we also require (9) to hold. This fixes for given a, 
£, B, the state vector up to a phase factor. In all consid- 
erations of flux densities only the quasi projectors enter 
and thus, such a phase factor in the states is irrelevant. 
The choice of phase factor is relevant for the calculation of 
the geometric phases via the line integrals (23) and (26). 
Then, we always make sure to choose a phase factor being 
differentiable along the path considered. 

Finally we note that for the case of no P violation, that 
is for 6 = 0, the numbering of the quasi projectors and the 
states \2Lj, F, F^, £, B) is done in a completely analogous 
way. 



B Relations for the geometric flux densities 

In this appendix we derive the representations (38), (39) 
and (40) for Ya^ab{K) and its derivatives, respectively. In 
the following we will omit the JiT-depcndence of all quan- 
tities for abbreviation. We now consider the expression 



Xr, 



■lb ■- 



i 1 
2 27ri 



Pn 



-Xi'^': 



Tr 



P 



f^ dKa ^'^ dKb 



{Q-Ep){C-E,f 



(B.l) 



According to the residue theorem the integral vanishes for 
/3 = 7 = a since a pole of third order aX C, = Ea gives a 
residual of zero 



^/^(C-i;^)^ 



27ri Rcs( 



(C - E^f 



■,C = E„) 



= . 
(B.2) 



Let Z? C (D be a simply connected set with Sa entirely 
inside D and E^, ^ D for all a ^ a; see Figure 2. Therefore, 
for 13 ^ a and 7 ^ a the integrand in (B.l) is analytic 
on £), and the integral in (B.l) vanishes due to Cauchy's 
integral theorem. The only two remaining cases 13 = a and 
7 7^ a as well as. j3 ^ a and 7 = a can be treated using 
again the residue theorem. Wc find easily 



\Y. 



2^ {Ea,- E^f 



-1 



Tr 



PaT7^P7T^ 



'dKa 'dKb 



2 f—^ {Ea — Ep) 



iy I 



rTr 



Pf 



-P 



dKa "dKb 



Tr 



T 



2£(^"-W 



dKa dKb 
dKb '^ dKa 



Tr 



5E 



1 



^i^S^--^^)' 



Tr 



_ d.^ ^ d.£ 
dKa dKb 



{a 4-)- b) 
(B.3) 



which is exactly Yaat, sec (37). Thus, we obtain the inte- 
gral representation (38) for Ya.ab 



^i£/^-[?c^ 



T. 



i 1 

227ri 



dCTr 




1 d^ 



C^- JLdKaiC,- J£)^ dKb\ 

(B.4) 



where we use the relation (15) for the quasi projectors in 
the last step. 
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In order to calculate the derivatives of Y^^ab we first we get 
derive some useful relations: 



d 



d 



[ic^^r'ic- 



dKa dKa 

^ [(C-^)-^](C-^) 



dKa 



<=> 



d 1 



1 



1 



dKa C - ^ C - ^ dKa C - ^ 



d 



-1 



d 



dKa dKa 
d 



[ic-^r'iC'^?] 



<^ 



^ 



[ic - ^r'] = ic - ^r 



dKa 



aifa ■ 2 27ri js^ 

1 52^ 1 

ir 



i< 



dyi 



(B.5) 



-Tr 
+ Tr 
-Tr 
-Tr 
+ Tr 



C - ^ dKadKb (C - ^)2 dK, 

1 a.^ 1 92.^ 



(C - ^)2 dKbC-^ dKadK^ 
1 d^ 1 d^ 1 



(9.^ 



C-^dKaC-^dKbiC- J£)^ dKc 

1 9^ 1 d^ 1 a.^ 



1 d^ 1 9^ 1 



1 a^ 1 9.-# 1 



C - ^ ^JsTb C - ^ c'/Ca (C - ^)2 dKc 



(B.8) 



9 



1 



1 



dKa (C - ^)2 (C - ^)2 dKaC-^ 

1 a^ 1 



c - ^ aA"a (c - ^y 



(B.6) 



Using the c 


yclicity of the trace this can be 


simplified to 


9if/-^-2 2.i£/^^{ 


Tr 


1 92^ 1 9^1 






C - .^ dKadKb (C - -^)2 ai^cJ 




+ Tr 


r 1 9^ 1 d.^ 1 a^i 


} 


U - ^ aif„ C--^dKh(C- ^)2 axJ 



- (6 f> c) , 



(B.9) 



With (B.5) and (B.6) we obtain from (B.4) 



which proves (39). With (15) we find 



5 _ i 1 

^"'"^ "2 2^ 



dKa 



Tr 



i< 



1 



d'^Ji 



1 



C-M. dKadKb (C - ^)2 dK^ 
1 dJ{ 1 d'^yi 



Tr 
Tr 
Tr 
Tr 



C^-JLdKbd- ^f dKadKb 
1 dJi 1 9^ 1 



C-^dKaC,~^dKb{C.- My dK, 
1 a.^ 1 dJ^ 1 



C - ^ 9^6 (C - ^)2 dKaQ-M dK, 
1 9^ 1 dJ^ 1 



Q-^dKbQ-M dKa (C 



aif. 






^{?£/^ 



(C-i^/3)-'(C-i?7)"'Tr 



d"^^ d.£ 



■ai^^aiffe ^9ife 



^ / dC(C-i?,^)^^(C-£^7)"'(C-i?.)' 



^ ,^ d.J( ^ d.J( ^ d./£ 



■[■ 



}-(&oc). (B.IO) 



~l The integrals in (B.IO) are easily evaluated using 

J • (i^-O Cauchy's theorems. 



With the short hand notations 



Using the cyclicity of the trace and performing partial « ._ p 

integrations of the second and fourth summand in (B.7) 



d^t 



dKa 



{ah) := 



d^ 



dKadKy 



(B.ll) 
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we obtain 
d 



dK 



y^M ^\{Y.^E^- ^^y^ Tr[a(a6)7c] 



With the short hand notation 
1 



7#a! 



C-^(5,B) 



5=0 



C-^^"'(^) 



(C.3) 



+ Y, -{E^ - Ep)-^ Tr[/3(a6)ac] 

- 2 ^ (£;„ - E^y^ TY[aaabac] 
-^^{Ea- E^y^ Tr[aa76ac] 
+ ^ (-B„ - Efiy^ Tr[^aa6ac] 

+ ^{E^- E^y\Ea - E^y TT[aa-/bac] 
+ ^ (^a - Ep)-\E^ - E,y TriPaabac] 

- ^ (-Ba - Ep)-\E^ - E^y TrilSajbac] 

/3,77^Q 

- H (-^" - E^y{E^ - Ep)-^ Tr[/3a76ac]} 

-(boc). (B.12) 

This can be simphfied to 

^Yo^M ^^{j2(Eo.' Epy Tr[a(a6)/3c - ac(3{ab)] 
+ 2. {Ea — Ep)^' Tr[— 2aaa6/3c + aafibac + I3aabac\ 



the expansion of the trace in (55) up to hnear order in the 
PV parameter 5 reads 



Tr 



1 



-D^ 



1 



:D, 



C - ^(£, B)—^ (c - M.{£, B)Y—'' 

= Tr[(^; + zS^^^z + 0{5'^))Dj 

X (^2 + z^S^^^z + z5^p^z^ + e'(5^))iZfc] 
= Tr[zZ2jZ^iZfe] + 5Tr[zi2j2^^Pv^iZfe 

+ zD_^z^-p^7?D_k + z^iZfe-z^pv^^iZj] + 0{8'^) . 

(C.4) 

Inserting (C.4) in (55) and performing a partial integra- 
tion we obtain 

+ (5Tr[z^pv zZJj-z^ZJfe - (j f^ fc)] 

(C.5) 

which proves (81). This derivation also holds for Z^ A. 
and I^ \J where i2iiiZfc are replaced by /^.,/i, and 
ZJj, M, J respectively. In this way we obtain from (56) and 
(57)"' 



/ ^ (£^a — Ep) (Ea — E^ 



r(B)PV^ 
a,jk 



i 1 



X Tr[aa/3bjc + /3aabjc — ja/Sbac — fiajbac] > 
-{b^c) (B.13) 

which proves (40). 

C Useful expressions for PV fluxes 

In this appendix we derive Iq, ,7^ (81) and give the anal- 
ogous expressions for I^ A. and X^ j^, . Then, we dis- 
cuss (70) and (71). 

Taking into account (2), (A. 5) and (A. 6) we first derive 
the expansion of (C,~ i^{.K))~'^ around J = 0. Analogously 
to (B.5) and (B.6) we find 



T).-riS.B) = 6--fJCT. 



-J^ 



C-^^°\£,By{C-^^°\£,B)) 



C,'^''\£,B)—^^ 

H,. - (.7 ^ k) 

(C.6) 



2tlk 



and 



lL',r''i£,B) = Si±£^dCTr 



1 



-^r 



-D, 



C-^(°H£,B)-'(C-^(°)(5,S)) 



C-^^'\e,B)"^'''' 

2l£k 



(C.7) 



d 1 



1 



1 



85 C-^ C-^ 85 C-^ 



(C.l) 



and 



1 



1 



8 1 ^ 

85 {C- ^)2 ~ (c-^)2~a5~c-^ 

1 8^ 1 



Now we discuss (70) and (71). From (69) we see that 
^ ^ '{£,B) and ^^ ' {—£, B) have the same set of eigen- 
values. We still have to check that our numbering scheme 
leads indeed to (70) and (71) for every a. We consider 
again only B ^ Q and vary £ starting from £ = 0: 



C-^ 85 {c,-^y 



(C.2) 



5(A) = X£ , 
A € [0, 1] 



(C.8) 
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For A =: (70) and (71) arc trivial. Increasing now A con- 
tinuously we encounter, due to B 7^ 0, no level crossings. 
Therefore, the identification of the eigenvalues and the 
quasi projectors corresponding to the same index a for 
(5(A), B) and (— 5(A),B) is always possible. This proves 
(70) and (71). 



where X3 = a/vA^ + (863(7/15)^. From the eigenvectors we 
calculate explicit representations of the projection opera- 
tors P?,^ for the 2S states and obtain 



D Detailed calculations of specific 
flux-density vector fields 

In this appendix we calculate the constants a" of (90). 
From (77) we find the P-violating part Ja (0, 3) of 

Ja (0, 3) to vanish. Therefore, neglecting terms of sec- 
ond order in the small PV parameter (5, we can calculate 
Ja (0, 3) setting (5 = 0. Then, the 2S states decouple 
from the 2P states, and we may restrict ourselves to the 
submatrix ^^^^°\q,3) oi ^^°'>{0,3) (2) with respect 
to the 2S states, see Tables 3 and 5 in Appendix E. The 
derivatives dBi^^^"^°Ho,3), i G {1,2,3}, of this subma- 
trix read in the basis \a,£ = 0,3 = 0) with a = 9, . . . , 12 
(see Table 2 in Appendix A) 



2S,(0) 



dBi 



2S,(0) 



dBo 



gfj-B 

2\/2 



2s/2 



dJ^ 



2S,(0) 



^ M, ' ' gpiB 



OBs 



'0 1 -1\ 
10 10 

10 1 
-10 1 0/ 

'0 -i i 

1 -i 
i i 
-i -i 

10 

1 

0-10 

,0100 



(D.l) 



Due to rotational invariance of J^(0, 3), see (90), we may 
specify 3 = B^e^ for the evaluation of a" in (91). This 
simplifies the calculation of the eigenvalues and of the 
right and left eigenvectors of ^^ ''°^(0, S). In this case 
we find 



^2S,(o)(o^^^g^) 



Xi 











Xi 
Xi 




(D.2) 



where xi = L+A/32-irs/2 and X2 = L-3A/32-irs/2. 
The eigenvalues of the matrix in (D.2) are 



Eg = Xi 
EiQ = L - 
Ell == Xi 



9^^B 



B. 



2 

^ X3 

32 16 

9l^B 



irs/2 



32 16 ' 



(D.3) 



P 



2S,(0) _ 



10 





,0000 



p 



2S,(0) 
10 



1 

2x^ 



Tp2S,(0) 



Tp2S,(0) 
^12 



/O 

X3+-4 



VO SBsgfiB 

/O 



10 

Vo 











SBsgi^B 

X3--4, 



1 

2x1 



'I 







Xs-^ 




Vo -SBsgfiB 




(D.4) 



Now, all ingredients for (37) are available, and a straight- 
forward calculation yields, with (42), (56), and (58), the 
P-conscrving flux-density vector field 



Jf(0,e = i?3e3) = 



■B3 63 

B3 es 

IK3P 





for a = 9 , 

for a ==11, (D.5) 

for a = 10, 12 . 



Rotational invariance of J^{0,3) then leads to (91). 



We now give the relations between the functions g", h", 
r = 1, . . . , 15, introduced in Section 4.3, and the functions 
gf, /if, i = 1,...,6, introduced for (117) and (118) in 
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Section 5.3: 
















g1 


= SoBo 


'g^ + i(2£| - £',)g- + i(2S| - 6^) g^ 




+ ^-1 


(D.6) 


~9^ 


= £„Bo [2B2£l 9ii + 2^2 qU , (D-7) 


5? 


= £oBo 


'-B2 gt + B2{£l - Ef) 5?! + B2 9^2 


, (D.8) 


r. 


= £oBo 


"-^ff6-^l.98-.9?o 


(D.9) 


~9t 


^ElBl 


'2 92 , (D.IO) 
53 , (D.ll) 
'h^ + i(2f| - £f) ht + \{2Bl - Bl) h^ 




+ \elBlK^ 


(D.12) 


^ 


= SoBo [2B2 K^ + 2B2BI /i?2] , (D.13) 


hs 


= foSo 


'-B2 K + B2{El - El) K, + B2EIBI K2 

(D 


5 

.14) 


K 


= foSo 


'-\BlK-K-Blh-, 


(D.15) 
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(D 
(D 
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E The matrix representations of ^q , D and /^ 

Tables 3, 4 and 5 show the mass matrix for zero external fields, ^ n, the electric dipole operator D_ and the magnetic 
dipole operator /i. for the n = 2 states of hydrogen. We give all these matrices in the basis of the pure 2S and 2P 
states, that is, the states for zero external fields and without the P-violating mixing. 
In Tables 4 and 5 we use the spherical unit vectors, which arc defined as 



1 , 
where Bi {i — 1, 2, 3) are the Cartesian unit vectors. For e±, the following relation holds: 



(E.l) 



(E.2) 



Table 3. The mass matrix ^q (3) for the n — 2 states of hydrogen. For the quantities A, L, A and 5\.2 see Table 1. The 
lifetimes rp,s are rs = T^ ^ = 0.1216s and tp = Fp^ = 1.596 x 10"^ s; see [37,38]. 





2P3/2,2,2 


2^3/2,2,1 


2^3/2,1,1 


2A/2,1,1 


2S'i/2,l,l 


2P3/2,2,2 


-hrp 














2^3/2,2,1 
2^3/2,1,1 

2A/2,1,1 
2Si/2,l,l 
















-krp 

A 




A 





+ ^52L 

L+^2 


192 ^/2 

t.-hrp 


192 v/2 




(Table 3a) 
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2P3 


/2,2,0 


2^3/2,1,0 


2Pi/2,l,0 


2Si/2,l,0 


2^1/2, 


0,0 


2Si/2,0,0 


2^3/2,2,0 





















2^3/2,1,0 







A- ^ 

^ 96 


A 













192 v^ 


2Pi/2,l,0 







A 


A i r 

96- 2^ P 


iSiL 

+^52L 










192%/2 


2S'i/2,l,0 










-iSiL 
-1S2L 


T , A 

^+32 










2Pi/2,0,0 
















A 
~32 - 


Kp 


iSiL 

+ |i52i 


251/2,0,0 
















-iSiL 
-|i<52L 







(Table 3b) 





2P3/2,2,-l 2P3/2,1,-1 2Pi/2,l,-l 2Si/2,l,-l 


2^3/2,2,-2 


2^3/2,2,-1 
2^3/2,1,-1 

2Pi/2,l,-l 
2Si/2,l,-l 


/I 1 -4 

i^i«' 

-^rp 

A -^ 

n 96 -^ n 








-^Fp 192v^ 
A A i^^ ^SiL 


^ 192V2 96 2^P +^S2L 

-2O2L -^±rs 


2^3/2,2,-2 


0000 





(Table 3c) 
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Table 4. The suitably normalised electric dipole operator 0/(ers) for the n = 2 states of hydrogen where rs is the Bohr 
radius for hydrogen. 





2P3/2,2,2 


2P3/2,2,1 


2^3/2,1,1 


2Pi/2,l,l 


2Si/2,l,l 


2P3/2,2,2 














-3e_ 


2^3/2,2,1 














75 ''o 


2^3/2,1,1 














-^eo 


2Pi/2,l,l 














-VSeo 


2Si/2,l,l 


3e+ 


75^" 


-\/ieo 


-\/3eo 






(Table 4a) 





2P3 


/2,2 





2P3 


/2,1,0 


2Pi/2,l,0 


2Si/2,l,0 


2Pi/2,0,0 


2Si/2,0,0 


2^3/2,2,1 

















-^- 








2^3/2,1,1 

















-\/l- 





-V6e- 


2Pi/2,l,l 

















-\/3e_ 





^/3e-. 


2Si/2,l,l 


V 


fie- 




^\ 


/§«- 


-VSe- 





VSe- 






(Table 4b) 





2^3/2,2,1 


2^3/2,1,1 


2Pi/2,l,l 


2Si/2,l,l 


2^3/2,2,0 











-\/i^+ 


2^3/2,1,0 











^/l-+ 


2A/2,1,0 











V3e+ 


2Si/2,l,0 


71^^+ 


\/l«+ 


V3e+ 





2Pi/2,0,0 











-\/3e+ 


25i/2,0,0 





V6e+ 


-^/3e+ 






(Table 4c) 





2^3/2,2,0 


2^3/2,1,0 


2Pi/2,l,0 


2Si/2,l,0 


2Pi/2,0,0 


25i/2,0,0 


2^3/2,2,0 











V&ea 








2^3/2,1,0 

















V6eo 


2Pi/2,l,0 

















-\/3e,) 


2Si/2,l,0 


VGso 











-\/3eo 





2Pi/2,0,0 











-\/3eo 








2Si/2,0,0 





\/6eo 


-\/3eo 












(Table 4d) 
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2^3/2,2,-1 


2^3/2,1,-1 


2Pi/2,l,-l 


2^1/2,1,-1 


2^3/2,2,0 











-\/l«- 


2^3/2,1,0 











-^/l-- 


2Pi/2,l,0 











-^/3e_ 


2Si/2,l,0 


^- 


-V^e- 


-VSe- 





2Pi/2,0,0 











-VSe- 


2Si/2,0,0 





VGe- 


-^/3e_ 






(Table 4e) 







2P3 


/2,2,0 


2P: 


/2, 1, 





2Pl/2,l, 





2Si/2,l,0 


2Pi/2,0, 





2Si/2,0,0 


2^3/2,2, 


-1 



















-75-+ 










2^3/2,1, 


-1 



















\/l-+ 







-V6e+ 


2A/2,1, 


-1 



















VSe+ 







VSe+ 


2Si/2,l, 


-1 


\ 


/ie+ 


V 


^e. 




^/3e+ 







^/3e+ 








(Table 4f) 





2^3/2,2,-1 


2^3/2,1,-1 


2A/2,1,-1 


2S'i/2,l,-l 


2^3/2,2, 


-2 


2^3/2,2,-1 











^-0 







2^3/2,1,-1 











\/l-o 







2A/2,1,-1 











VSeo 







2Si/2,l,-l 


^-0 


yl^o 


VSeo 





3e_ 




2^3/2,2,-2 











-3e+ 






(Table 4g) 
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Table 5. The suitably normalised magnetic dipole operator /x//is for the n = 2 states of hydrogen, where /is = eh/{2m^) is 
the Bohr magneton and g — 2.002319304(76) is the Lande factor of the electron [63]. 





2P3/2,2,2 


2P3/2,2,1 


2^3/2,1,1 


2A/2,1,1 


2Si/2,l,l 


2^3/2,2,2 


-^eo 


v^(9+2) „ 
4 ^-^ 


v^(9+2) „ 
4^/3 


-^- 





2^3/2,2,1 
2^3/2,1,1 

2Pi/2,l,l 
2Si/2,l,l 











f7f-° 







-feo 



(Table 5a) 



2^3/2,2,0 



2^3/2,1,0 



2Pi/2,l,0 2Si/2,l,0 2Pi/2,0,0 



2Si/2,0,0 



2^3/2,2,1 
2^3/2,1,1 

2Pi/2,l,l 
25'i/2,l,l 



73(9 + 2) , 


4 ' 


9 + 2 
12 ^- 


^(9-1) f. 



9+2p 

5(9 + 2) „ 
12 '^ 

^(9-1) . 



^(9-1) , 
6 







^(9-1) , 
3 



(Table 5b) 



2^3/2,2,1 



2P3/2,1,1 2Pi/2,l,l 2Si/2,l,l 



2P3/2,2,0 

2P3/2,1,0 

2Pi/2,l,0 
2Si/2,l,0 
2Pi/2,0,0 
2Si/2,0,0 



v^(9+2) „ 
4 ^+ 



-S±2e4 

4^3^^ 



9^g 






12 '^+ 



5(9 + 2) 


12 '^ 

%/2(9-l) 


6 





e+ 



v^(9-i) , 

6 ' 

^(9-1) , 
6 '^ 

9-4„ 

-■Ve+ 



9-4„ 
6-6+ 









§«+ 



2 + 



(Table 5c) 



2P3/2,2,0 



2P3/2,1,0 2Pi/2,l,0 25'i/2,l,0 2Pi/2,0,0 2Si/2,0,0 



2P3/2,2,0 

2P3/2,1,0 

2Pi/2,l,0 
2Si/2,l,0 
2Pi/2,0,0 
2Si/2,0,0 





9+2 



eo 



V2(9-l) 
3 







eo 



9+2 
6 





eo 






V^(9-l) 
3 





eo 



-^^^eo 




9-4 



eo 









■^^^^eo 

9-4 







eo 



(Table 5d) 
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2^3/2,2,-1 



2P3/2,1,-1 2Pi/2,f,-f 



25'i/2,l,-l 



2^3/2,2,0 

2^3/2,1,0 

2Pi/2,l,0 
2Si/2,l,0 
2Pi/2,0,0 
2Si/2,0,0 



Va(9+2) , 

4 

9 + 2 „ 
4^3*^- 







a+lp 

12 *=- 

5(g + 2) 
12 



6 








V2(3-l) , 
3 ' 





6 " 

v^(9-l) , 
6 ' 

6 ^- 









(Table 5e) 



2^3/2,2,0 



2P3/2,1,Q 2Pi/2,l,Q 2Si/2,l,0 2Pi/2,0,0 25i/2,0,0 



2P3/2,2,-l 



2P 



3/2 



1, 



2Pi/2,l,-l 
2Si/2,l,-l 



%/3(9+2) 


4 ' 


9+2 
12 '^ 


^(9-1) 



6 ^+ 





4v^'^+ 

5(9+2) 
12 '^+ 

6 ^+ 





V^(9-l) 



e+ 



9-4„ 

-Ve+ 







2 + 





V^(9-l) 
3 

9-4 



e+ 



6 e + 





(Table 5f) 





2P3/2,2,-l 2P3/2,1,-1 2Pi/2,l,-l 2S'i/2,l,-l 


2P3/2,2,-2 


2P3/2,2,-l 
2P3/2,1,-1 

2Pi/2,l,-l 
2Si/2,l,-l 


^-0 -f±|eo -^eo 
-f$i«o %2)eo -l^eo 

-^«o -§71'^° -"^^^ 
feo 


V^(9+2) „ 
4 '^- 

V5(9+2) „ 
4v^ ^- 




2P3/2,2,-2 


.2(1+2),^ vi^,^ _^,^ 


2±^eo 



(Table 5g) 



T. Gasenzer et al.: Metastable states of hydrogen and their geometric phases and flux densities 



29 



References 

1. M.V. Berry, Proc. R. Soc. Lond. A392, 45 (1984) 

2. B. Simon, Phys. Rev. Lett. 51, 2167 (1983) 

3. A. Shapere, F. Wilczek, eds.. Geometric Phases in Physics, 
Vol. 5 of Advanced Series m Mathematical Physics (World 
Scientific, Singapore, 1989) 

4. M. Nakahara, ed.. Geometry, Topology and Physics, Grad- 
uate student series in physics (Adam Hilger, Bristol and 
New York, 1990) 

5. J. Garrison, E. Wright, Physics Letters A 128, 177 (1988) 

6. S. Massar, Phys. Rev. A 54, 4770 (1996) 

7. F. Keck, H.J. Korsch, S. Mossmann, Journal of Physics A: 
Mathematical and General 36, 2125 (2003) 

8. M.V. Berry, Czechoslovak Journal of Physics 54, 1039 
(2004) 

9. W. Heiss, Czechoslovak Journal of Physics 54, 1091 (2004) 

10. A.I. Nesterov, F.A. de la Cruz, Journal of Physics A: Math- 
ematical and Theoretical 41, 485304 (2008) 

11. T. Bergmann, T. Gasenzer, O. Nachtinann, Eur. Phys. J. 
D 45, 197 (2007) 

12. T. Bergmann, T. Gasenzer, O. Nachtmann, Eur. Phys. J. 
D 45, 211 (2007) 

13. M. DeKieviet, D. Dubbers, C. Schmidt, D. Scholz, 
U. Spinola, Phys. Rev. Lett. 75, 1919 (1995) 

14. T. Bergmann, M. DeKieviet, T. Gasenzer, O. Nachtmann, 
M.L Trappe, Eur. Phys. J. D 54, 551 (2009) 

15. M. DeKieviet, T. Gasenzer, O. Nachtmann, M.L Trappe, 
Hyperfine Int. 200, 35 (2011) 

16. LB. Khriplovich, Parity Nonconservation in Atomic Phe- 
nomena (Gordon and Breach, Philadelphia, 1991) 

17. M. Bouchiat, C. Bouchiat, Rep. Prog. Phys. 60, 1351 
(1997) 

18. M.A. Bouchiat, [physics.atom-ph] arXiv:1111.2172vl 
(2011) 

19. M.A. Bouchiat, J. Guena, L. Hunter, L. Pottier, Phys. 
Lett. B 117, 358 (1982) 

20. C.S. Wood, S.C. Bennett, D. Cho, B.P. Masterson, J.L. 
Roberts, C.E. Tanner, C.E. Wieman, Science 275, 1759 
(1997) 

21. S.C. Bennett, C.E. Wieman, Phys. Rev. Lett. 82, 2484 
(1999) 

22. M.J.D. Macpherson, K.P. Zetie, R.B. Warrington, D.N. 
Stacey, J. P. Hoare, Phys. Rev. Lett. 67, 2784 (1991) 

23. N.H. Edwards, S.J. Phipp, P.E.G. Baird, S. Nakayama, 
Phys. Rev. Lett. 74, 2654 (1995) 

24. P.A. Vetter, D.M. Meekhof, P.K. Majumder, S.K. Lamore- 
aux, E.N. Fortson, Phys. Rev. Lett. 74, 2658 (1995) 

25. D.M. Meekhof, P. Vetter, P.K. Majumder, S.K. Lamore- 
aux, E.N. Fortson, Phys. Rev. Lett. 71, 3442 (1993) 

26. K. Tsigutkin, D. Dounas-Frazer, A. Family, J.E. Stalnaker, 
V.V. Yashchuk, D. Budker, Phys. Rev. Lett. 103, 071601 
(2009) 

27. K. Nakamura et al. (Particle Data Group), J. Phys. G 37, 
075021 (2010) 

28. L.W. Wansbeek, B.K. Sahoo, R.G.E. Timmermans, 
K. Jungmann, B.P. Das, D. Mukherjee, Phys. Rev. A 78, 
050501 (2008) 

29. O.O. Versolato et al.. Physics Letters A 375, 3130 (2011) 

30. A. Schafer, G. Soff, P. Indelicato, B. Miiller, W. Greiner, 
Phys. Rev. A 40, 7362 (1989) 

31. L.N. Labzowsky, A.V. Nefiodov, G. Plunien, G. Soff, 
R. Marrus, D. Liesen, Phys. Rev. A 63, 054105 (2001) 



32. V.M. Shabaev, A.V. Volotka, C. Kozhuharov, G. Plunien, 
T. Stohlker, Phys. Rev. A 81, 052102 (2010) 

33. F. Ferro, A. Artemyev, T. Stohlker, A. Surzhykov, Phys. 
Rev. A 81, 062503 (2010) 

34. F. Ferro, A. Surzhykov, T. Stohlker, Phys. Rev. A 83, 
052518 (2011) 

35. R.W. Dunford, R.J. Holt, J. Phys. G: Nucl. Part. Phys. 
34, 2099 (2007) 

36. R.W. Dunford, R.J. Holt, Hyperfine Int. 200, 45 (2011) 

37. J. Sapirstein, K. Pachucki, K.T. Cheng, Phys. Rev. A 69, 
022113 (2004) 

38. L.N. Labzowsky, A.V. Shonin, D.A. Solovyev, J. Phys. B: 
At. Mol. Opt. Phys. 38, 265 (2005) 

39. V.F. Weisskopf, E.P. Wigner, Z. Phys. 63, 54 (1930) 

40. V.F. Weisskopf, E.P. Wigner, Z. Phys. 65, 18 (1930) 

41. L. Mandel, E. Wolf, Optical coherence and quantum optics 
(Cambridge University Press, Cambridge, 1995) 

42. S.M. Barnett, P.M. Radmore, Methods in Theoretical 
Quantum Optics (Oxford University Press Inc., New York, 
1997) 

43. P.L. Knight, L. Allen, Phys. Lett. A 38, 99 (1972) 

44. Y.K. Wang, I.C. Khoo, Opt. Comm. 11, 323 (1974) 

45. W. Bernreuther, O. Nachtmann, Z. Phys. A 309, 197 
(1983) 

46. O. Nachtmann, Elementary particle physics, concepts and 
phenomena (Springer, Berlin, 1990) 

47. G.W. Botz, D. Bruss, O. Nachtmann, Annals of Physics 
240, 107 (1995) 

48. R. Jacob, R.G. Sachs, Phys. Rev. 121, 350 (1961) 

49. R.G, Sachs, Annals of Physics 22, 239 (1963) 

50. H. Flanders, Differential Forms, Vol. 11 of Mathematics 
in Science and Engineering (Academic Press, New York, 
1963) 

51. D. Bruss, T. Gasenzer, O. Nachtmann, Eur. Phys. J. direct 
D2, 1 (1999) 

52. C.E. Loving, P.G.H. Sandars, J. Phys. B: At. Mol. Phys. 
10, 2755 (1977) 

53. M. DeKieviet, private communication 

54. S.F. Pate, D.W. McKee, V. Papavassiliou, Phys. Rev. C 
78, 015207 (2008) 

55. D. de Florian, R. Sassot, M. Stratmann, W. Vogelsang, 
Phys. Rev. Lett. 101, 072001 (2008) 

56. D. de Florian, R. Sassot, M. Stratmann, W. Vogelsang, 
Phys. Rev. D 80, 034030 (2009) 

57. M. Alekseev et al. (COMPASS), Physics Letters B 693, 
227 (2010) 

58. E. Leader, A.V. Sidorov, D.B. Stamenov, [hep-ph] 
arXiv:1103.5979v2 (2011) 

59. U. Jentschura, S. Kotochigova, E. LeBigot, P. Mohr, 
B. Taylor, The Energy Levels of Hydrogen and Deuterium, 
National Institute of Standards and Technology, Gaithers- 
burg, MD (2005) 

60. S.G. Karshenboim, Phys. Rep. 422, 1 (2005) 

61. J. Erler, M.J. Ramsey-Musolf, Phys. Rev. D 72, 073003 
(2005) 

62. E.U. Condon, G.H. Shortley, The Theory of Atomic Spec- 
tra (Cambridge University Press, 1963) 

63. P.J. Mohr, B.N. Taylor, Rev. Mod. Phys. 77, 1 (2005) 



